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Abstract 

This paper is devoted to the study of couplings of the Lebesgue measure and the Pois- 
son point process. We prove existence and uniqueness of an optimal (i.e. 'asymptotically 
optimal' and 'translation invariant') coupling whenever the asymptotic mean transportation 
cost is finite. Moreover, we give precise conditions for the latter which demonstrate a sharp 
^ threshold aX d — 2. The cost will be defined in terms of an arbitrary increasing function of 

the distance. 

The coupling will be realized by means of a transport map ('allocation map') which 
assigns to each Poisson point a set ('cell') of Lebesgue measure 1. In the case of quadratic 
costs, all these cells will be convex polyhedra. 

P_i 1 Introduction and Statement of Main Results 

+^ a) Given a translation invariant point process fi* : uj fi^ = X]5es(aj) ^(0 ' with 

rH unit intensity, we consider the set H of all couplings q* of the Lebesgue measure £ and the point 

I " I process - i.e. the set of measure-valued random variables io s.t. for a.e. to the measure q^ 

on X M"^ is a coupling of £ and fi"^ - and we ask for a minimizer of the asymptotic mean cost 
functional 



i?(|x - y\)dq'{x,y) 



> 1 
I£;i Coo(g') := liniinf— — E 

00 

Here Bn := [0,2")'^ C W^. The scale i? : M+ ^ M+ will always be some strictly increasing, 
CN continuous function with ■d{Q) = and lim 7?(r) = oo. 

o 

A coupling w i— )■ of the Lebesgue measure and the point process is called optimal if it 
minimizes the asymptotic mean cost functional and if it is translation invariant in the sense 
"'"j that its distribution is invariant under push forwards of the measures dq^{x,y) on x by 

translations (x, y) i— )• (x + y + z), z G 'L'^. Our main result states 



Theorem 1.1. // the asymptotic mean transportation cost 

^{\x - y\)dq*{x,y) 



liminf inf — -K 

n-s>oo g«en £(i3„) 



'x_B„ 



is finite then there exists a unique optimal coupling of the Lebesgue measure and the point process 

b) The unique optimal coupling q^ can be represented as {Id,T'^)^,£, for some map T'^ : M"^ — t- 
supp(//'^) C M*^ measurably only dependent on the sigma algebra generated by the point process. 
In other words, defines a (deterministic) fair allocation rule. Its inverse map assigns to each 
point ^ of the point process ('center') a set ('cell') of Lebesgue measure ^'^{S,) G N. If the point 
process is simple then all these cells have volume 1. In the case of quadratic cost, i.e. i?(r) = r^, 
the cells will be convex polyhedra. The transport map will be given as = V(p^ for some 
convex function 99"^ : M'^ — t- M and induces a Laguerre tessellation (see |LZ08j ). 
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In the case '!?(r) = r the transportation map induces a Johnson-Mehl diagram (see |Aur91] ). 
For the many results on and apphcations of these tesseUations see the references in [LZ08j and 
|Aur9 1] . In the hght of these results one might interpret the optimal coupling as a generalized 
tessellation. 

c) As a particular corollary to Theorem 1.1 we conclude that Cqo = infq«gn ^oo{q*) and that the 
infimum is always attained, more precisely, it is attained by a translation invariant coupling q*. 
For translation invariant couplings q' the mean cost functional j^IE /jgd^yi ^(1^ ~ vl) dq'{x, y) 
however, is independent of A C M'^. Hence, 



^inf E 



d{\x-y\)dq*{x,y) 

'x[0,l)'* 



where liinv now denotes the set of all translation invariant couplings of the Lebesgue measure 
and the point process. 

Moreover, for translation invariant couplings, the mean cost of transportation E [&{\x — T'{x)\)] 
of a Lebesgue point x to the center of its cell is independent of x G W^. Hence, 

Coo=infE[^?(|0-r'(0)|)] 

where the infimum is taken over all translation invariant maps T : M'^ x — )■ with ^2, = 
for a.e. w. And again: the infimum is attained by a unique such T. 

d) Analogous results will be obtained in the more general case of optimal 'semicouplings' 
between the Lebesgue measure and point processes of 'subunit' intensity. 

We develop the theory of optimal semicouplings as a concept of independent interest. Optimal 
semicouplings are solutions of a twofold optimization problem: the optimal choice of a density 
p < 1 of the first marginal and subsequently the optimal choice of a coupling between p^i 
and /i2- This twofold optimization problem can also be interpreted as a transport problem with 
free boundary values. 

Given a point process of subunit intensity and finite mean transportation cost we prove that 
there exists a unique optimal semicoupling between the Lebesgue measure and the point process. 
It can be represented on X as before as q^ = {Id,T'^)^:2, in terms of a transport map 
T'^ : M"^ — )■ supp[^'^] U {9} where S now denotes an isolated point ('cemetery') added to W^. 

e) In any case, we prove that the unique transport map can be obtained as the limit of a 
suitable sequence of transport maps which solve the optimal transportation problem between 
the Lebesgue measure and the point process restricted to bounded sets. 

More precisely, for z ^TL^ and 7 S F := ({0, 1}'^)^ consider the 'doubling sequence' of cubes 

n 

B^{z,j) = z-^2'=-V. + [0,2")'='. 

k=l 

Note that the cube Bn{z,^) is one of the subcubes obtained by subdividing Bn+i{z,j) into 2^^ 
cubes of half edge length. Let Tz^n{-,^,j) '■ — ^ supp[^'^] U {S} be the transport map for the 
unique optimal semicoupling between £ and 1b„(^,^) • /x"^, that is, for the optimal transport of 
an optimal 'submeasure' p'^ ■ £, to the point process restricted to the cube Bn{z,^). 

Theorem 1.2. For every z £ Z'^ and every bounded Borel set M <ZW^ 

\\m{£.®¥(^v){{{x,u),-i) e M xVLxT : Tzn{x,uj,-f) / T{x,uj)}) = 

where v denotes the Bernoulli measure on F. 
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f) If /u* is a Poisson point process with intensity /3 < 1 we have rather sharp estimates for the 
asymptotic mean transportation cost to be finite. 

Theorem 1.3. (i) Assume d > 3 (and f3 < 1) or /3 < 1 (and d > 1). Then there exists a 
constant < k. < oo s.t. 

. logT?(r) ^ ^ log??(r) 

lim mf -j — < K Coo < oo =^ hm sup — < k. 

r— 5>oo r"- r— >co 

(ii) Assume d <2 and (3 = 1. Then for any concave i? : [l,oo) — )• M dominating •& 

'dir) 'Qir) 

-dr<oo Coo < oo =^ hmsup , = 0. 



1 I ' I r— >oo ' ' 

The first imphcation in assertion (ii) is new. Assertion (i) in the case /? = 1 is due to Hohoyd 
and Peres |HP05j . based on a fundamental result of Talagrand |Tal94j . The second implication 
in assertion (i) in the case /? < 1 was proven by Hoffman, Holroyd and Peres |HHP06] . The 
second implication in assertion (ii) is due to |HL01| . 

Now let us consider the particular case of transportation cost, i.e. i?(r) = r^. 

Corollary 1.4. (i) For all d € N, all (3 < 1 and p £ (0,cxd) the asymptotic mean L^- 
transportation cost Coo is finite if and only if 

_ _ ( oo, for d > 3 or f3 < 1; 
^^^'^ { i, for d<2 and (3 = 1. 

(ii) If (3 = 1 then for allp G (0, oo) there exist constants < k < k' < oo s.t. for all d > 2(p A 1) 

k-dP/^ < Coo < k'-dP/^. 



g) The study of fair allocations for point processes is an important and hot topic of current 
research, see e.g. |HP051 ITim08| IHPPS09] and references therein. A landmark contribution was 
the construction of the stable marriage between Lebesgue measure and an ergodic translation 
invariant simple point process |HHP06] . One of the challenges is to produce allocations with 
fast decay of the distance of a typical point in a cell to its center or of the diameter of the cell. 
The gravitational allocation |CPPR06l ICPPRj in d > 3 was the first allocation with exponential 
decay. Moreover, all the cells are connected and contain their center. However, the decay was 
not yet as good as the decay of a random allocation constructed in |HP05| . 

On the other hand, during the last decade the theory of optimal transportation (see e.g. 
|RR98j . |Vil03j ) has attracted lot of interest and has produced an enormous amount of deep 
results, striking applications and stimulating new developments, among others in PDEs (e.g. 
|Bre91j . |Ott01| . [XGS08 ]). evolution semigroups (e.g. |OV00| . |AZ09) . |OS09j ) and geometry 
(e.g. (StuOeilEtuOebj, [D7n9] . (Vn09], |Ohtn9j V Ajtai, Komlos and Tusnady as weh as Tala- 
grand and others studied the problem of matchings and allocation of independently distributed 
points in the unit cube in terms of transportation cost ( |AKT84] . [Tal94j and references therein). 

h) In all the transportation problems considered in the afore mentioned contributions, however, 
the marginals have finite total mass. Our paper seems to be the first to prove existence and 
uniqueness of a solution to an optimal transportation problems for which the total transportation 
cost is infinite. 

More precisely, the main contributions of the current paper are: 

• We present a concept of 'optimality' for (semi-) couplings between the Lebesgue measure 
and a point process. Even in the particular case of semicouplings between the Lebesgue 
measure and a finite counting measure, this concept is new. 
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Figure 1: Optimal semicoupling of Lebesgue and 25 points in the cube with cost function 
c{x, y) = \x — and (from left to right) p=l, 2, 4 respectively. 

• We prove existence and uniqueness of an optimal semicoupling whenever there exists a 
semicoupling with finite asymptotic mean transportation cost. 

• We prove that for a.e. doubling sequence of boxes {Bn{z,'y))^^j^ the sequence of optimal 
semicouplings qn,z,'y between the Lebesgue measure and the point process restricted to the 
box Bn{z,^) will converge. More precisely, the sequence g*^^,^ will converge as n — oo 
towards the unique optimal semicoupling q' between the Lebesgue measure and the point 
process. 

• We prove that the asymptotic mean transportation cost for the Poisson point process in 
d < 2 is finite for L^-costs with p < d/2 or, more generally, e.g. for 'd{r) = r'^l'^ ■ [y^^^^Y 
with a > 1. 

2 Set-up and Basic Concepts 

£ will always denote the Lebesgue measure on M"^. 
2.1 Couplings and Semicouplings 

For each Polish space X (i.e. complete separable metric space) the set of measures on X - 
equipped with its Borel a- field - will be denoted by Given any ordered pair of Polish 

spaces X, y and measures A S A4(X),/i G M.{Y) we say that a measure g G A^(X x y) is a 
semicoupling of A and briefly q £ ns(A,/i), iff the (first and second, resp.) marginals satisfy 

(7ri)*g < A, {iT2)*q = /i, 

that is, iff q{A x Y) < X{A) and q{X x B) = ^{B) for aU Borel sets A C X,B C Y. The 
semicoupling q is called coupling, briefly q G n(A, /u), iff in addition 

(vri)*g = A. 

Existence of a coupling requires that the measures A and fj, have the same total mass. If the total 
masses of A and // are finite and equal then the 'renormalized' product measure q = j^x)^ ^ 
is always a coupling of A and ^. 

If A and ^ are S-finite, i.e. A = Yl^=i^n-: = Yl^=il^n with finite measures A„ G M{X), 
G A4{Y) - which is the case for all Radon measures - and if both of them have infinite total 
mass then there always exists a S-finite coupling of them. (Indeed, then the A„ and /i„ can be 
chosen to have unit mass and q = X]n('^" does the job.) 
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2.2 Point Processes 

Throughout this paper, /x* will denote a translation invariant point process of subunit intensity, 
modeled on some probability space (r2,2l, P). For convenience, wc will assume that S7 is a 
compact metric space and 21 its completed Borel field. These technical assumptions are only 
made to simplify the presentation. 

Recall that a point process is a measurable map //* : ^ A^(M'^), u: ^ ijl^ with values in the 
subset Mccmnt^'^) of locally finite counting measures on W^. It is a particular example of a 
random measure, characterized by the fact that jJL^^^A) G Nq for P-a.e. oj and every bounded 
Borel set A(ZW^. It can always be written as 

with some countable set H(a;) C M'^ without accumulation points and with numbers k{^) G N. 
The point process is called simple iff k{^) = 1 for all ^ G '^(^) and a.e. lo or, in other words, iff 
IJ,{{x}) G {0, 1} for every a; G M*^ and a.e. lo. 

The point process /j,* will be called translation invariant iff the distribution of fx* is invariant 
under push forwards by translations : x x + z of M'', that is, iff 

for each z G M'^. We say that n* has subunit intensity iff E [/i*(yl)] < 2>{A) for all Borel sets 
A c M°'. If "=" holds instead of "<" we say that /x* has unit intensity. A translation invariant 
point process has subunit (or unit) intensity if and only if its intensity 

is < 1 (or = 1, resp.). 

Given a point process /Lt*, the measure d(/x*P)(?/, w) := diJ,^{y) dF{oj) onW^ x CI will be called 
universal measure of the random measure /x*. 

The most important example of a translation invariant simple point process is the Poisson point 
process or Poisson random measure with intensity fi It is characterized by 

• for each Borel set A C M*^ of finite volume the random variable ui i->- ii'^{A) is Poisson 
distributed with parameter j3 ■ S,{A) and 

• for disjoint sets Ai, . . . A^ dW^ the random variables ij!^{Ai), . . . , iJi!^{Ak) are independent. 

There are some instances in which we need additional assumptions on fi* (e.g. ergodicity, unit 
intensity). In each of these cases we will clearly point out the specific assumptions we make. 



2.3 Couplings of the Lebesgue Measure and the Point Process 

A (scmi-)coupling of the Lebesgue measure £ G M.{W^) and the point process /x* : f2 — >■ M.{W^) 
is a measurable map q* : Q. M.{W^ x W^) s.t. for P-a.e. u 

(f is a (semi-) coupling of Z and [i^ . 

We say that a measure Q G M.{^^ x M'^ x 0) is an universal (semi-) coupling of the Lebesgue 
measure and the point process iff dQ{x, y, cu) is a (semi-)coupling of the Lebesgue measure dii(x) 
and of the universal measure d{fi*F){y, uj). 

Disintegration of an universal (semi-)coupling w.r.t. the third marginal yields a measurable map 
q* : ^ A^(M'^ X W^) which is a (semi-)coupling of the Lebesgue measure £ and the point 
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process ji* . Conversely, given any (semi-) coupling q* of the Lebesgue measure 2, and the point 
process then 

dQ{x,y,u) := dg"(x, y)dP(a;) 
defines an universal (semi-)coupling. 

According to this one-to-one correspondence between q* — (semi-)coupling of £ and //* — and 
Q = q*W — (semi-)coupling of £ and /i*P — we will freely switch between them. In many cases, 
the specification 'universal' for (semi-) couplings of 2, and /x*P will be suppressed. And quite 
often, we will simply speak of (semi-) couplings of S, and fi*. 

2.4 Cost Functionals 

Throughout this paper, will be a strictly increasing, continuous function from R+ to M+ with 
t?(0) = and lim t?(r) = oo. Given a scale function -& as above we define the cost function 

c(x,y) = -d{\x- y\) 

on X W^, the cost functional 

Cost(g) = / c{x,y)dq{x,y) 

on M. {W^ X R^) and the mean cost functional 

£ost((5) = / c{x,y) dQ{x,y,uj) 

on M. {W^ xW^ xQ,). We have the following basic result on existence and uniqueness of optimal 
semicouplings the proof of which is deferred to the Appendix. 

Theorem 2.1. (i) For each bounded Borel set A dW^ there exists a unique semicoupling Qa 
of 2, and (1a/x*)P which minimizes the mean cost functional Cost(.). 

(ii) The measure Qa can he disintegrated as dQA{x,y,uj) := dq^(x,y) dF{uj) where for ¥ -a. e. uj 
the measure q'^ is the unique minimizer of the cost functional Cost(.) among the semicouplings 
of 2 and IaA*'^- 

(Hi) £05t(QA) = J^Cost{q'X)dP{oo). 

For a bounded Borel set A C M'' , the transportation cost on A is given by the random variable 
Ca ■ ^ ^ [0, oo] as 

Ca{u)) := Cost(g^) = inf{Cost(g'^) : q'^ semicoupling of £ and 1a m'^}- 

Lemma 2.2. (i) If Ai, . . . , A n are disjoint then Vo; G O 

n 

C n ^ (W) > ^Ca,(w) 

i=l «=1 

(ii) If Ai and A2 are translates of each other, then Cai and C^j are identically distributed. 

(Hi) IfAi,...,A n are disjoint and /i*(Ai), . . . ,/x*(^n) are independent, then the random vari- 
ables CAi,i = 1, . . . , n, are independent. 
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Figure 2: Concept of exhausting sequences: start with a smah cube and repeatedly double its 
edge lengths to exhaust space (cost function c{x,y) = \x — yp). 

Proof. Property (ii) and (iii) follow directly from the respective properties of the point process 
and the invariance of the Lebesgue measure under translations. The intuitive argument for (i) 
is, that minimizing the costs on [J - Ai is more restrictive than doing it separately on each of the 
Ai. The more detailed argument is the following. Given any semicoupling g'^ of £ and 
then for each i the measure := iRdxA,?'^ is a semicoupling of £ and lyi^/i'^. Choosing q'^ as 
the minimizer of C n (w) yields 

U 

i=l 

Cu^A,(^) = Cost(g-) = 5^Cost(gr) >J2^aM- 

i i 

□ 

2.5 Convergence along Standard Exhaustions 

For n € No := N U {0} and z G Z*^ define the cube or box Bn{z) of generation n with basepoint 
z by 

Br^iz) = z + [0,2^y. 
For z = simply put 5„ = B„(0). More generally, for 7 G F := ({0, l}"^)^ put 

n 

Bn{z,j) = z-^2^-Sfc + [0,2")'^. 

k=l 

Starting with the unit box Bq{z,j) = 2; + [0, 1)'^, for any random vector 7 G F the sequence 
{Bn{z,j))neNo can be constructed iteratively as follows: Given the box Bn{z,^) attach 2*^ — 1 
copies of it - depending on the random variable 7n+i = (7n+i! • • • 5 7n+i) with values in {0, 1}'^ - 
either on the right (if 'Jn+i = 0) or on the left (if 7^+1 = 1), either on the backside (if j'^+i ~ 0) 
of on the front (if j^+i — either on the top (if Jn+i — 0) or on the bottom (if 7^^^ = 1), etc. 
The sequence {Bn{z,'y))ne^o fixed z and 7 is increasing and for i/-almost every 7 G F it 
increases to M'^. Each of the boxes Bn{z,"f) contains the point z. 

Put 

c„ := 2""'" • E [Cb„(z,^)] . 

Note that translation invariance implies that the right hand side does not depend on z G Z"^ and 
7 G F. 
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Corollary 2.3. (i) The sequence {Cn)n^N() 'is non- decreasing. The limit 



lim c 



sup c„ 



exists in (0, oo]. 

(a) Assume that /i* is ergodic. Then, we have for all z G 1/^ , for all ^ and for ¥ -almost 
every u; G 17; 



Iiminf2-"'^C 



,7)('^) 



(Hi) Coo ^ iiif £00(9) where Tig denotes the set of semicouplings of 2 and fi* . 

Proof, (i) is an immediate consequence of the previous lemma. For (ii) fix an arbitrary nested 
sequence of boxes {Bn)n generated by a standard exhaustion. Then we have by superadditivity 
Vw G 17 for all n, fc G N 



where Bn are disjoint copies of -B„ such that Uj=i = Bn+k- In the limit of /c — )• 00 we get by 
ergodicity for P-a.e. 00 

-kdr 



for each n G N and thus 



liminf2-'='^CB,(w) > E 

fc— ^-oo 



liminf 2-'='^Cb,(w) > Cc 



n—nd/~ 



On the other hand, Fatou's lemma implies 



E 



< liminf E 



-nd 



Iiminf2-"°'CB 

_ n— >oo 

Both inequalities together imply the assertion. 

For (iii) take any semicoupling q* of £ and /i*P. Then, we have for any n 

Taking the limit yields 

^00(9) = liminf 2"'^"e:ost(lj^dxS„xfi9*) > limc„ = Coo- 

n— s>oo n 

□ 

Corollary 2.4. Cqo only depends on the scale and on the distribution of n' , - not on the 

choice of the realization of on a particular probability space (0,21, P). 

Proof. It is sufficient to show, that Cn just depends on the distribution of fi*. For a given set 
of points ^'(0;) in Bn there is a unique semicoupling q'^ of £ and lB„fJ''^ minimizing Cost (see 
Proposition 6.3). Hence, just depends on ^"(0;). However, the distribution of the points in 
Bn, ^i^jj), just depends on the distribution of □ 

Remark 2.5. All the previous definitions and results did not require that has subunit 
intensity. However, one easily verifies that 



/3 > 1 



00 



where /3 := E [;U*([0, 1)'^)] denotes the intensity of the translation invariant point process. 
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Figure 3: Semicoupling of Lebesgue and 25 points in the cube with c(x, y) = \x — y\ where each 
point gets mass 1/9, 1/3, 1 respectively. 

Remark 2.6. The problem of finding an optimal semicoupling between £ and a Poisson point 
process fi* of intensity /3 < 1 is equivalent to the problem of finding an optimal semicoupling 
between £ and (3 ■ fT where fi' is a Poisson point process of unit intensity. 
Indeed, given /3 E (0, 1) and a semicoupling q* of £ and a Poisson point process ^i* of intensity 
13. Put T : X ^ f3~^/'^ X on as well as on x W^. Then [i^ := t^,^^ is a Poisson point process 
with intensity 1 and 

r ■■= P ■ r*g" 

is a semicoupling of £ and f^-fi^. Conversely, given any Poisson point process fi^ of unit intensity 
any any semicoupling q'^ of £ and (3 ■ jl^ then q'^ := ^ ■ {T~'^)^,q'^ is a semicoupling of £ and 
:= {T~^)^:jl'^, the latter being a Poisson point process of intensity (3. In both cases, q is 
translation invariant if and only if q is translation invariant. 

The asymptotic mean transportation cost for q* measured with scale ?? will coincide with the 
asymptotic mean transportation cost for q* measured with scale '&p{r) := (3 ■ 'd{f3~^^'^ r): 

E [ ^{\x - y\) dq' =E [ - y\) dq\ 

3 Uniqueness 

Throughout this section we fix a translation invariant point process A4(K'^) of subunit 

intensity and with finite asymptotic mean transportation cost Cqo- 

Proposition 3.1. Given a semicoupling q'^ of £ and fi^ for fixed w G then the following 
properties are equivalent: 

(i) For each bounded Borel set A C M'^, the measure l^dxyl'?'^ ^-^ unique optimal semicou- 
pling of the measures '■= £(.) — q^{.,ZA) and Ia/^*^- 

(Note that the measure Xj^ is the part of £ which - under the coupling q^ - is not transported 
to li^^/x*^. If q^ is a coupling then the definition simplifies to X^{-) = q^{.,A), see Figure 

w 

(a) The support of q^ is c-cyclically monotone, more precisely, 

N N 
'^c{xi,yi) < ^c(xi,yi+i) 

1=1 i=l 

for any N £ N and any choice of points . . . , {xN,yN) in supp(g'^) with the con- 

vention yN+i = yi- 
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Figure 4: The left picture is a semicoupling of Lebesgue and 36 points with cost function 
c{x,y) = \x — y\^. In the right picture the five points within the small cube can choose new 
partners from everything that is white (corresponding to the measure A';^). If the semicoupling 
on the left hand side is locally optimal, then the points in the small cube on the right hand side 
will choose exactly the partners they have in the left picture. 



(in) There exists a set C M'^ and a c- cyclically monotone map : — )• such that 

^- = (/(i,r-),(u.£). (3.1) 

Recall that, by definition, a map T is c-cyclically monotone iff the closure of its graph 
{{x,T{x)) : X £ A'^} is a c-cyclically monotone set. 



Proof. The implications (Hi) =^ (ii) =^ (i) follow from Lemma 6.1 



(i) =^ (Hi): Fix an exhaustion (S^)n of M*^ by boxes, say B'^ = [—2" ^,2" ^)^. For each n G 
N, let p';^ be the density of the measure X'^ := A*^, on M*^. This is the part of Lebesgue measure 
from which the points inside of B'^ might choose their 'partners'. Obviously, < < Pn+i < 1 
and / 1 a.e. on M'=' as n ^ oo as A^^(.) = £(.) - q'^{., CM'^) = £. 

Assuming (i), according to Proposition 6.3 (or, more precisely, a canonical extension of it for 
semicouplings of /)£ and a) there exists a set A':^ and a c-cyclically monotone map : — )■ M'^ 
such that 

dq^{x,y) = dbT^i^^){y)\At{x)pl{x)dZ{x) on M'^ x S^. 

Since the left hand side is independent of n, obviously p^{x) G {0,1} a.s., that is, without 
restriction p%^{x) = 1. This in turn implies A^ C A^j^^ (up to sets of measure 0) and 



J-n+l — J^n on 



This trivially yields the existence of 

:= lim on A"^ := hm <, 

defining a a c-cyclically monotone map A"^ — )• M'^ with the property that 

d4^{x,y) = d6T^i^^){y)lA'-{x)d2.{x). 

□ 

Remark 3.2. In the sequel, any transport map : A^ — ?• as above will be extended to a 
map T'^ : M'^ — 7- M'^ U {9} by putting T'^(x) := S for all x £W^\ A^ where S denotes an isolated 
point added to ('point at infinity', 'cemetery'). Then (3.1) simplifies to 



q'^ = {Id, T'^)^ £ onW X W. (3.2) 



10 



Definition 3.3. . 

\> A semicoupling Q = q*F of £, and fi* is called locally optimal iff some (hence every) of 
the properties of the previous proposition are satisfied for P-a.e. a; G Jl. 

> A semicoupling Q = q*F of £, and fi* is called asymptotically optimal iff 

liminf 2~"'''(2:ost(l]B;dxB' Q) = ^oo 
for some exhaustion {B'^)n o/M'^ by boxes B'^ = Bn{z,j). 

> A semicoupling Q = q*F of Z and fi* is called translation invariant iff for each z G Z*^ the 
measure Q is invariant under the translation 

{x,y,uj) t-^ {x + z,y + z,uj) 

ofW^ xR^xn. 

\> A semicoupling Q = q'F of & and fj,* is called optimal iff it is translation invariant and 
asymptotically optimal. 

The very same definition of 'optimality' applies to couplings of & and fj,*. 

Remark 3.4. (i) Asymptotic optimality is not sufficient for uniqueness and it does not imply 
local optimality: Given any asymptotically optimal semicoupling q* and a bounded Borel set 
^ C M'^ of positive volume, choose an arbitrary coupling q^^ of the measures q'^{.,A) and Ia^"^ 
— which are the marginals of := Ir^xA^'^- If M'^(^) ^ 2 (which happens with positive 
probability) then one can always achieve that q'^ is a non-optimal coupling and that it is different 
from q'^. Put 

q ■=Q +<iA-(iA- 

Then q' is an asymptotically optimal semicoupling of £ and ji* . It is not locally optimal and it 
docs not coincide with q* . 

(ii) Local optimality does not imply asymptotic optimality and it is not sufficient for unique- 
ness: For instance in the case p = 2, given any coupling q* of £ and /x* and z G M*^ \ {0} 
then 

dcf{x,y) := dq'^ix + z,y) 

defines another locally optimal coupling of £ and /x* . At most one of them can be asymptotically 
optimal. 

(iii) Note that local optimality — in contrast to asymptotic optimality and translation invari- 
ance — is not preserved under convex combinations. We do not claim that local optimality and 
asymptotic optimality imply uniqueness. 

Given 7, 77 G M.{W^) with ^{W^) > r){W^) we define the transportation cost by 

Cost(7,r?) := inf{Cost(g) : g G ^s(7,^?)}• 
Similarly, given measure valued random variables 7* , ry* : O — >• (M"^) and a bounded Borel set 
A C M'' we define the mean transportation cost by 

e:ost(7*,r?*) := inf {e:ost(g*P) : q"^ G Usi^f"^ ,ri'') for a.e. u} . 
Given a coupling Q = q*F of £ and /i*P we define the efficiency of the coupling Q on the set A 

by 

^A[Q) ■= — ZT, 7^- 

£0St(lRdxA<V) 

It is a number in (0, 1]. The coupling Q is said to be efficient on A iff eA{Q) = 1- Otherwise, it 
is inefficient on A. 
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Lemma 3.5. (i) Q is locally optimal if and only if ca{Q) = 1 for all bounded Borel sets A C M'^. 
(a) ca{Q) = 1 for some A cW^ implies CA'iQ) = 1 for all A' C A. 

Proof, (i) Let A be given and w € be fixed. Tlien l^d^^Q^ the optimal semicoupling of the 
measures A';^ and IaI^-^ if and only if 



Co5t{l^,^AQn = Cost (A^ , lAli'^) . 
On the other hand, CAiQ) = 1 is equivalent to 

E [Cost(lRdxA(?')] = E [Cost (A^ , 



(3.3) 



The latter, in turn, is equivalent to (3.3) for P-a.e. u; G fi. 

(ii) If the transport q restricted to x A is optimal then also each of its sub-transports. □ 

Theorem 3.6. Every optimal semicoupling of £, and ^'P is locally optimal. 

Proof. Assume we are give a coupling Q of £ and fx*¥ which is translation invariant and not 
locally optimal. According to the previous lemma, the latter implies that there exist n G N and 
zq S Z'^ such that the coupling Q is not efficient on the box i?„(zo), i-e. 

By translation invariance this implies ^b„{z){Q) = rj < 1 for all z £ Z'^. Hence, for each z € Z*^ 
there exists a measure- valued random variable 9^,1(2) such that Qb„{z) is a semicoupling 

of A"^ ^^-j and 1b,^(z)IJ''^ and more efficient than ^^-^ := l]Rtix_B„(2) ' l'^^ i-^- such that 



E 



Cost(^^, 



< ri-E 



Cost(g|;, 



n{z)) 



Put 



ze{2" 



Then q* is a semicoupling of £ and fi* and for all z € (2" Z)'' 



E 



CoSt(l]udxB„(2)Q'* 



< ?7-E 



Costfl. 



<B„{z)Q 



Translation invariance of q* - together with uniqueness of cost minimizers on bounded sets - 
implies translation invariance of q* under the group (2"Z"'). In other words, Q = q'F IS a 
translation invariant semicoupling of £ and /i*P which satisfies 

(2:0St(ljjdxR < V ^05t{l^dy,B„(z)Q) 



for all z G (2"Z)'^. Additivity of the mean cost functional (iosi{.) implies 

e:OSt(lRdxB„+, 0) < ?7-£0St(lKdxB„+fe<3) 

for all G No and therefore, due to Corollary |2.3[ iii), finally 

Coo < liminf Cost(l]f5dxB, 0) < • liminf £05t(liRdxB. Q) 

fc— 5>oo fc— >oo 

with r] < 1. This proves that Q is not asymptotically optimal. 

Theorem 3.7. There exists at most one optimal semicoupling of £ and //'P. 



□ 
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Proof. Assume we are given two optimal semicouplings and q'. Then also q* := ^Q* + -^Q* 
is an optimal semicoupling. Hence, by the previous theorem all three couplings - q' and 
q* - are locally optimal. Thus, for a.e. uj by the results of Proposition 3.1 there exist maps 
T^,T^,T'^ and sets Af, A^,A^ such that 

d6rp^i^^){y)lA'-{x)d2.{x) = dq''{x,y) 

^d6T^^^){y)lA-^{x) + ^dST^(^^){y)lA^{x)^ d2{x) 

This, however, implies T^{x) = T2{x) for a.e. x E n A2 and, moreover, A^ = A2. Thus 
q^ = q^. □ 

Remark 3.8. Note that we only used translation invariance under the action of Z"^. However, 
the minimizer is translation invariance under the action of M'^. For the uniqueness it would also 
have been sufficient to require translation invariance under the action of kU^ for some /c E N. 

Theorem 3.9. (%) If fi* has unit intensity then every optimal semicoupling of 2 and fi* is indeed 
a coupling of them. 

(a) Conversely, if an optimal coupling exists then fi* must have unit intensity. 

Proof, (i) Let Q be an optimal semicoupling. For n G N put Bn{z) = z + [0,2")'^ and consider 
the saturation ak '■= 2~^'^Q{Bk{z) x Bk{z) x Q) < 1. Note, that ak is independent of 2: G Z"^. 
Hence, we have ak < a^+i. Indeed, Bk+i{z) is the disjoint union of 2^^ cubes Bk{yj) for suitable 
yj. Therefore, 

2d 

ak+1 > 2-'^^2-^'='Q(5fc(yj) x Bk{yj) x n) = a^. 
i=i 

Thus, the limit Ooo := hnifc-^-oo Ofc exists and we have aoo £ (0, 1]. 

Since n' has unit intensity and since Q is a semicoupling we have Q(M'^ x Bk x il) = 2^'^. Let us 
first assume that aoo < 1 and choose r = [(1 + ^(1 — aoo))^^^ — l]/2. Then for all G N mass of 
a total amount of at least (1 — aoo)2^'^ has to be transported from CB^ into B^. The volume of 
the (r2'')-neighborhood of the box Bj. is less than ^(1 — aoo)2^'^. Hence, mass of total amount 
of at least 5(1 — aoo)2^'^ has to be transported at least the distance r2''. Thus, we can estimate 
the costs per unit from below by 

2-"" [ c{x,y)dQ{x,y,u) > 1(1 - a^)^(r2'=). 

The right hand side diverges as k tends to infinity which contradicts the finiteness of the costs 
per unit. Thus, we have Ooo = 1- Furthermore, for all k there is a u G Bk{0) such that 

ak = 2~''^Q{Bk{0) X Bk{0) X n) = 2"'=^ ^ Q{Bo{v) x BkiO) x n) 

< Q{Bq{u) X Bk{0) xn) < Q{Bq{u) xW^ X O). 

However, by translation invariance the quantity Q{Bq {u)xR'^x 0) is independent of u. Moreover, 
it is bounded above by 1 as Q is a semicoupling. Hence, we have for all v G M'^: 

1 = limsupoA: < Q{Bo{v) x R'^ x Q) < I. 

k—^oo 

Therefore, Q is actually a coupling of the Lebesgue measure and the point process. 

(ii) Assume that Q is an optimal coupling and that /3 < 1. Then a similar argumentation as 
above yields that for each box Bk, Lebesgue measure of total mass > (1 — /?) • 2^'^ has to be 
transported from the interior of Bk to the exterior. As k tends to 00, the cost of these transports 
explode. □ 
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Corollary 3.10. In the case ??(r) = r^, given an optimal coupling q* of 2, and a point process 
/X* of unit intensity then for a.e. £ Q there exists a convex function (^"^ : M'^ — )• M (unique up 
to additive constants) such that 

q^ = {Id,V^n.^- 

In particular, a 'fair allocation rule' is given by the monotone map = Vif^. 
Moreover, for a.e. u and any center £^ G "^iy^) '■= supp(/i'^), the associated cell 

is a convex polyhedron of volume G N. // the point process is simple then all these cells 

have volume 1. 



Proof. By Proposition 3.1 we know that T"^ = lim„_s.oo Tni where is an optimal transporta- 
tion map from some set to B^. From the classical theory (see |Bre9H[GM96j ) we know that, 
= V(^J^ for some convex function y?^. More precisely, 

</.^(x)= max {x" -\x-i\^ /2 + bi.) 

for some constants b^. Moreover, we know that = on for any € N. Fix any 

Co G Then, there is n e N such that e Then, (r;^+fc)~HCo) = {TnY^Ho) for any 

A; € N. Furthermore, 

Ti;i{x) = io ^ -|x-eolV2 + %, >-k-eiV2 + &€ veeHHni?;, e/Co- 

For fixed this equation describes two halfspaces separated by a hyperplane (defined by 

equality in the equation above). The set S^{£,q) is then given as the intersection of all these 
halfspaces defined by and ^ G H(a;) n Hence, it is a convex polytope. □ 

4 Construction of Optimal Semicouplings 

Again we fix a translation invariant point process /i* : — )■ A^(M"') of subunit intensity and 
with finite asymptotic mean transportation cost Cqo- 

4.1 Second Randomization and Annealed Limits 

The crucial step in our construction of an optimal coupling of Lebesgue measure and the point 
process will be the introduction of a second randomization, — besides the first randomness 
modeled on the probability space (r2,2t,P) which describes the random choice w i— t- /x*^ of a 
realization of the point process. The second randomization describes the random choice 7 1— )• 
{Bnjz, 'y))„ fzn o f a-ii increasing sequence of boxes containing a given starting point z G Z*^ (see 
also section 2.5). It is modeled on the Bernoulli scheme (r,5S(r), u) with F = ({0, l}*^)^, 5S(r) 
its Borel cr-field and v the uniform distribution on F = ({0, l}"^)^ (or, more precisely, the infinite 
product of the uniform distribution on {0, l}'^). 

For each z G G F and A; G N, recall that QBk{z,'y) denotes the minimizer of (tosi among the 

Translation invariance of 



2.1 



semicouplings of £ and (1^^.(2,^) /U*)P as constructed in Theorem 
this minimizer implies that d(5Bj.(z',^)(x, y, w) = dQB^.{z,'y)ix + z — z' ,y + z — z' ,uj) for all 
z,z' G Z'^. Put 



dQ^^{x,y,uj) := J dQB^(z,'y)ix,y,uj)diy{j) 



and dQ';{x,y,uj) := lBo(2)(y)dQ^(x, y, w). 

The measure defines a semicoupling between the Lebesgue measure and the point process 
restricted to be box Bq{z). It is a deterministic, fractional allocation in the following sense: 
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• it is a deterministic function of /i^ and does not depend on any additional randomness 
(coming e.g. from dv^'j)) 

• the measure transported into a given point of the point process has density < 1. 

The last fact of course implies that the semicoupling is not optimal. The first fact implies 
that all the objects derived from in the sequel - like QJ^ and Q°° - are also deterministic. 

Lemma 4.1. (i) For each k eN and z eZ'^ 

J c{x,y)dQ'^{x,y,uj) < Coo- 
R<^xBo{z)xn 

(a) The family (<52)feeN of probability measures onW^ xW^ xQ, is relatively compact in the 
weak topology. 

(Hi) There exist probability measures Q"^ and a subsequence {ki)i^fi such that for all z G Z*^; 

Qx' — > Of weakly as Z ->■ oo. 

Proof, (i) Let us fix ^; G Z'' and start with the important observation: For given n G N the 
initial box Bq{z) has each possible 'relative position within Bn{z,j) ' with equal probability. 
Hence, together with translation invariance of QBk{z,'y) (which in turn follows from that of P) 
we obtain 

j c{x,y)dQ'^{x,y,uj) = J J c{x,y)dQB^^(z,'Y){x,y,oj)du{j) 



ixBo{z)xa r mdxBo{z)xn 



veBk{z 

2—kd 



J c{x,y)dQBi^{z){x,y,u) 



R<^xBoiv)xn 



J c{x,y)dQBk{z){x,y,uj) 



'R<i'xBk{z)xQ. 

(ii) In order to prove tightness of ((5^)fceN) let 

:= {yeM.^: inf \x - y\ < m} 

xeBoiz) 

denote the closed m- neighbor hood of the unit box based at z. Then 

Q^(Cir^ X So(^) X 0) < [ c(x,y)dQ^(x,y,u;) 

v[m) jRdy^Bo{z)xa 

1 

Since — > oo as m — )■ oo this proves tightness of the family (Qz)fceN on M"^ x M*^ x fi. (Recall 
that 0, was assumed to be compact from the very beginning.) 

(iii) Tightness yields the existence of and of a converging subsequence for each z. A stan- 
dard argument ('diagonal sequence') then gives convergence for all z G Z'^ along a common 
subsequence. □ 
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Lemma 4.2. (i) For each r > there exist numbers £fe(r) with £fe(r) — >■ as k oo such that 
for all z, z' e 'L^ and allk eN 

Jr Jr 7 

for any Borel set A C'R''' xR'^ x ft. 

(a) For all zi, . . . ,Zm & Z*^, allk eN and all Borel sets A C M'' 

m / \ 

1=1 V 1=1 / 

Proof (i) Firstly, note that for each z,z' eZ'^,k eN,^ eT: 

z' G Bk{z, 7) ^ 3y : 7) = Bk{z', 7') 

and in this case 

uiW : Bk{z',j') = Bk{z,j)}) = 2-^-'. 

Moreover, 

: z' ^ Bk{z,^)]) < eki\z-z'\) 

for some efe(r) with £/;(r) — >■ as A; — >■ 00 for each r > 0. It imphes that for each pair z, z' G 'L'^ 
and each A; G N 

z/({7 G r : 37' : ^^(^,7) = i?fc(/,7')}) > 1 - ^k{\z - A)- 
Therefore, for each Borel set A C M'' x M*^ x ^^ 

/ QB^{z',-i){A)dv{a) < / QBk{z,^)iA)di^(.7) + £k{\z - z'\) ■su■pQB^(z',^){A)■ 
J^ Jr 7 

(ii) According to the previous part (i), for each Borel set ^ C M'^ 

m 

^Q^.(AxM<^x 0) 

i=l 

m „ 

= J2 QB,{zi,j){AxBo{zi)xn)du{-f) 
i=l -^^ 

< y] / Qsfc(^i,7)(^ X Bo{zi) X Q)di/(7) + ek{\zi - Zi\) ■ supQs^(^,,^)(^ x BQ{zi) x 
1^1 L-^r 7er 

m 

1=1 

□ 

Theorem 4.3. The measure Q°° : — X^zeZ'* optimal semicoupling of 2, and ji' . 
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Proof, (i) Second/third marginal: For each Borel set A C M x 



(vr2 3)*Q'^(M'^ X A) = yQf{M.'^xA)= Y lim Q^J {M.^ x A) 

^ — ' ^ — ' Z— >-oo 
= 2^/1™/ lBo(^)'^'3B,,(.,7)'^r(7) 

= y lim {fi'F){A n (5o(2) X Q)) 
(ii) First marginal: Let an arbitrary bounded open set A C M'^ be given and let (zi)igN be an 



enumeration of Z . According to the previous Lemma 4.2, for any m S N and any k £N 



i=l \ i=l 

Letting first k tend to oo yields 



1=1 

Then with m — ?■ oo we obtain 
which proves that (7ri)*Q°° < £. 



^Q^(v4 X M'^ X fi) < £{A). 
1 

Q^{A xR"^ xn)< £.{A) 



(iii) Optimality: By construction, Q°° is translation invariant. Due to its translation invariance, 
the asymptotic cost is given by 



/ c{x,y)dQ'^{x,y,u;) = / c{x , y) dQf {x , y , uj) 

c{x,y)dQ^{x,y,uj) < Coo- 



'x_Bo(0)xQ 



Here the final inequality is due to Lemma 4.1, property (i) (which remains true in the limit 



k = oo), and the last equality comes from the fact that 

/ c{x,y)dQ^^{x,y,ui) = 

for all z ^ u and for all G N (which also remains true in the limit k = oo). □ 
Corollary 4.4. (i) For k — )■ oo, the sequence of measures := ^ Q^, /c G N, converges 
vaguely to the unique optimal semicoupling Q°° . 

(ii) For each z & the sequence ((5^)fceN converges vaguely to the unique optimal semicoupling 

Proof, (i) A slight extension of the previous Lemma 4.1(iii) + Theorem 4.3 yields that each 
subsequence (Q'^")n of the above sequence {Q^)k will have a sub-subsequence converging vaguely 
to an optimal coupling of 2, and fi*. Since the optimal coupling is unique, all these limit points 
coincide. Hence, the whole sequence (Q'^)fc converges to this limit point (see e.g. |Dud02j . Prop. 
9.3.1). 
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(ii) Lemma 4.2 (i) implies that for z, z', u ^T,^ and every measurable A C 

< Efcdz - z'l) • sup Qb^(v){A ^ ^o(w) X 0)) 

< efc(|z-z'|)^0 

as A: — )• oo. Hence, for each f €^ Cc{ Q) and each z' G R 



x,y,uj)dQ^ 



0. 



That is, \jfdQ^-j f dQ';,\ ^ as A: ^ oo. 



□ 



Corollary 4.5. ftave Coo = inf Coo(^*) where Hg denotes the set of all semicouplings q* 
of £, and jj,* . In particular, it holds that 



inf lim inf — r 



E 



/ c{x,y)dq'{x,y) 



1 



lim inf inf 

rn-oo q'sHs £(i3„ 



-E 



c{x,y) dq'{x,y) 



Proof. The optimal coupling Q constructed in the previous Theorem has mean asymptotic trans- 
portation cost bounded above by Cqo- Thus, we have inf Coo(<?*) < Cqo- Together with Lemma 

2.3, this yields the claim. □ 



4.2 Quenched Limits 

According to chapter 3, the unique optimal semicoupling between d£{x) and d^^{y) dF{Lo) can 
be represented on x M'^ x 17 as 

dQ'^(x,y,io) = d5Ticc,u.){y) d£{x) d¥{oj) 

by means of a measurable map 



X n 



u{g}. 



defined uniquely almost everywhere. Similarly, for each z ^ Z'^ and A; G N there exists a 
measurable map 

r^,fc : M'^ X 17 X r ^ U {9} 
such that for each 7 S F the measure 

dQBk{z,^){x,y,i^) = d5T^ i^{x,uj,^){y) d£.{x)dF{u) 

on M'^ X M'^ X 17 is the unique optimal semicoupling of d£,{x) and 15^.(2,7) (y) dji^iy) dF{uj). 

Proposition 4.6. For every z ^ Z'^ 

Tz^k{x,io,^) — )• T{x,uj) as A: — )• 00 in £, 01^ -measure. 

The claim basically relies on the following lemma which is a slight modification (and extension) 
of a result in |Amb03j . 

Lemma 4.7. Let X,Y be locally compact separable spaces, 9 a Radon measure on X and p a 
metric on Y compatible with the topology. 

(i) For alln £N let Tn,T : X ^ Y be Borel measurable maps. Put dQn{x, y) := d5xj^(x)iy)d9{x) 
and dQ{x,y) := d5x[x){y)d9{x). Then, 



Tn ^ T in measure on X 



Qn Q vaguely in M-{X x Y). 
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(n) More generally, let T and Q as before whereas 

dQn{x,y) := / d6T„Ux'){y)de'{x')d9{x) 
Jx' 

for some probability space (X',2l', 9') and suitable measurable maps r„ : X x X' — > y. Then 

Qn ^ Q vaguely in M.{X x Y) ^> Tn{x,x') — >■ T{x) in measure on X x X' . 

Proof, (i) Assume r„ — > T in ^-measure. Then also / o [Id, Tn) f o (Id, T) in ^-measure for 
any f ^ Cc{X xY). Therefore, by the dominated convergence theorem we have 

j f{x,y)dQn = I f{x,Tn{x))de ^ I f{x,T{x))de = J f{x,y)dQ. 

This proves the vague convergence of Qn towards Q. 

For the opposite direction, fix K C X compact and £ > 0. By Lusin's theorem there is a compact 
set if C K such that T\k is continuous and 0{K\K) < e. Put rj :R+ ^ R+,t t-^ 1 A \t\ /e. The 
function 

(f>{x,y) = iK{x)v{p{y,T{x))) 

is lower semicontinuous, nonnegative and compactly supported. Hence, there exist (pi G Cc{X x 
Y) with 4>i \ (j). By assumption, we have for each I 



j (t>{x,y)dQn{x,y) < j (j)i{x,y)dQn{x,y) J (j)i{x,y)dQ{x,y). 



Moreover, 



J (t>l{x,y)dQ{x,y) J (t){x,y)dQ{x,y) = 0. 
Therefore, limn^oo / (f'ix, y)dQn{x, y) = 0. In other words, 

lim / lK{x)r]{p{Tn{x),T{x)))de{x) = Q. 



This implies lim„_).oo 0{{x G K : p{Tn{x),T{x)) > e}) = and then in turn 

lim e{{x G K : p{Tn{x), T{x)) > 2e}) = 0. 

n— >oo 

(ii) Given any compact K C X and any e > 0, choose d) as before. Then vague convergence 
again implies lim„^.oo / <P{x,y)dQn{x,y) = 0. This, in other words, now reads as 

lim / / lK{x)r]{p{Tn{x,x'),T{x)))de'{x')de{x) = Q. 
Jx Jx' 

Therefore, 

lmi^{9 9') (^i^^{x,x') e K X X' : p{Tn{x,x'),T{x)) > 2e^^ =0. 
This is the claim. □ 

Proof of the Proposition. Fix z £ Z'^ and recall that 

^ Q°° vaguely on R'^ x M.'^ 

where 

dQ^ix,y,u) = dST^x^^){y)d£{x)dF{uj) 

and 

dQ'l{x,y,u) = j^dQB^{z,^){x,y,u:)du{-i) = J^d6T^^^^x,oo,j)iy) d&ix) dF{uj) di^ij) 

with transport maps T : x O ^ M'^ U {9} and T^^k : x x F ^ M<^ U {9} as above. 
Apply assertion (ii) of the previous lemma with X := x U,X' = F, y = M*^ U {9} and 

e = £(g)F,0' = u. □ 



19 



Actually, this convergence result can significantly be improved. 
Theorem 4.8. For every z ^ and every bounded Borel set M C M*^ 

lim(£(g)P(g)z/)({(x,a;,7) gMkOxT: T^^k{x,u,-i) ^ T{x,u)}) = 0. 

fe— >-oo ' 

Proof. Let M as above and £ > be given. Finiteness of the asymptotic mean transportation 
cost implies that there exists a bounded set M' C such that 

(£®P) ({(ar,a;) G M X : T{x,uj) ^ M']) < e. 

Given the bounded set M' there exists 5 > Q such that the probability to find two distinct 
particles of the point process at distance < 5, at least one of them within M', is less than e, i.e. 

P {[u : 3{y, J/') G M' X M<^ : Q<\y-y'\<5, ^x''{{y}) > 0, /^'^({y'}) > o}) < e. 

On the other hand, Proposition 4.6 states that with high probability the maps T and T^^^ have 
distance less than 5. More precisely, for each (5 > there exists such that for all k > ko 

(£(8)P(g)z/)({(a;,a;,7) G M x O x T : \T^^k{x,0J,j) - T{x,oj)\ > S}) < e. 

Since all the maps T and T^^k take values in the support of the point process (plus the point S) 
it follows that 

(£(g)P(g)z/)({(x,a;,7) gMxOxT: T^^k{x,oj,-i) ^ T{x,u)]) <^e 

for all k > ko- □ 

Corollary 4.9. There exists a subsequence {ki)i such that 

Tz^ki{x,(^,"y) — ^ T{x,u)) as / — >■ oo 

for almost every x G M.'^, G 11, 7 G F and every z G Z'^. Indeed, the sequence {Tz^ki)l is finally 
stationary. That is, there exists a random variable : M*^ x O x F ^ N such that almost surely 

Tz,ki{x,uj,j) = T{x,uj) for all Z > /^(x, 7). 

Corollary 4.10. There is a measurable map T : M{W^) M{R'^ x W^) s.t. q"^ := T(/x'^) de- 
notes the unique optimal semicoupling between S, and fi'^ . In particular the optimal semicoupling 
is a factor coupling. 

Proof. The map oj ^ q^ is measurable with respect to the sigma algebra generated by /i* . Thus 
there is a measurable map T such that q* = □ 

5 Estimates for the Asymptotic Mean Transportation Cost of a 
Poisson Process 

Throughout this section, n* will be a Poisson point process of intensity /3 < 1. The asymptotic 
mean transportation cost for /x* will be denoted by 

Coo = Coo{'&,d,P) 

or, if 'd{r) = r^, by Coo{p,d,f3). We will present sufficient as well as necessary conditions for 
finiteness of Coo- These criteria will be quite sharp. Moreover, in the case of iT'-cost, we also 

present explicit sharp estimates for Cqo- 

To begin with, let us summarize some elementary monotonicity properties of Coo(^?, d, 
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Lemma 5.1. (i) 'd < ■& implies Coo{'&,d,f3) < Coo{'&,d, (3). 

More generally, limsupj._i.oo < ^jo and Coo{'&,d,f3) < oo imply Coo{'d,d,l3) < oo. 

(ii) If = 99 o -i? for some convex increasing ip : M-|_ — t- M-|- then ip (^l3~^Coo{i}, d, 13)^ < 
^-h^{^,d,(3). 

(Hi) 13 < [3 implies Coo{^ ■,d^j3) < Coo('!?, /?)• 

Proof, (i) is obvious, (ii) If q denotes the optimal semicoupling for ■!? then Jensen's inequality 
implies 

r'coo(^,d,/3) = / if{^{\x-y\))dq{x,y) 

J]Rdx[0,l)'* 

> iflr^Ef i^ilx - y\) dq{x,y)] >ip{p-hoo{^,d,p)). 

(iii) Given a realization Jtf^ of a Poisson point process with intensity /3. Delete each point 
^ G supp[/I'^] with probability 1 — 13/(3, independently of each other. Then the remaining point 
process /i"^ is a Poisson point process with intensity /3. Hence, each semicoupling q^ between £ 
and Ji^ leads to a semicoupling q^ between £ and /x"^ with less or equal transportation cost: the 
centers which survive are coupled with the same cells as before.) □ 



5.1 Lower Estimates 



Theorem 5.2 ( ^HLOlj ). Assume (3=1 and d < 2. Then for all translation invariant couplings 
of Lebesgue and Poisson 



E 



'x[0,l)'* 



\x -y 



2'ix,y) 



00. 



Theorem 5.3. For all (3 < 1 and d > 1 there exists a constant k' = n'{d,(3) such that for all 
translation invariant semicouplings of Lebesgue and Poisson 



E 



'x[0,l)'' 



exp ( k'|x — yl'^ ) dq*{x,y) 



00. 



The result is well-known in the case (3 = 1. In this case, it is based on a lower bound for the event 
"no Poisson particle in the cube [— r, r)'^" and on a lower estimate for the cost of transporting 
the Lebesgue measure in [— r/2,r/2)'^ to some distribution on M'^ \ [— r, r)*^: 



Coo > exp 



(2r) 



..(-). 2- 



Hence, Cqo — ^ 00 as r — )• 00 if t?(r) = exp(K;' r'^) with k' > 2^°^. 

However, this argument breaks down in the case (3 < 1. We will present a different argument 
which works for all /3 < 1. 

Proof. Consider the event "more than (3r)'^ Poisson particles in the box [— r/2, r/2)'^" or, for- 
mally, 

n{r) = {f^'{[-r/2,r/2)'^)>{3rf}. 

Note that E/x* ([-r/2, r/2)'^) = /3r^ with p < 1. For CO G ^{r), the cost of a semicoupling 
between £ and l[_r./2,r/2)d/^'^ bounded from below by 



t?(r/2) ■r'^ 
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(since r'^ Poisson points - or more - must be transported at least a distance r/2). The large 
deviation result formulated in the next lemma allows to estimate 

for any k < //^(S'^) and suitable r„ — t- oo. Hence, if '!9(r) > exp(K'r'^) with k' > 2'^ ■ //3(3'^) then 

Coo > ^i^{rn)) ■ i9(r/2) • > exp((K'2"'^ - ky^) ■ r'^ ^ oo 
as r — )■ oo. □ 
Lemma 5.4. Given any nested sequence of boxes Bn{z,j) C M*^ 



n— ^-oo Z 



2nd 



lj'{Bn{z,^))>t 



with Ip{t) = tlog{t//3) -t + 13. 



Proof. For a fixed sequence Bn{z,^), n G N, consider the sequence of random variables Zn{.) = 
lj,'{Bn{z,^)). For each n G N 

with Xi = fj,*{Bo{i)). The Xi are iid Poisson random variables with mean (3. Hence, Cramer's 
Theorem states that for all t > /3 



liminf 7-^ log I 



1 

2nd 



Zn > t 



> m 



with 



/^(t) = sup [tx - log/i(x)] = tlog(t//3) - t + p. 



□ 



5.2 Upper Estimates for Concave Cost 

In this section we treat the case of a concave scale function -i?. In particular this implies that 
the cost function c{x,y) = i^^x — y\) defines a metric on M'^. The results of this section will be 
mainly of interest in the case d <2; in particular, they will prove assertion (ii) of Theorem 1.3. 
It suffices to consider the case /3 = 1. Similar to the early work of Ajtai, Komlos and Tusnady 
[AKT84j . our approach will be based on iterated transports between cuboids of doubled edge 
length. 
We put 

0{r) := / i?(s)ds and e{r) := sup , . 

Jo s>r S'^l^ 



5.2.1 Modified Cost 

In order to prove the finiteness of the asymptotic mean transportation cost, we will estimate the 
cost of a semicoupling between £ and 1^^* from above in terms of the cost of another, related 
coupling. 

Given two measure valued random variables z^* , z/* : — )■ 7W(M"') with z^^(M"') = z/^(M'^) for a.e. 
cj G we define their transportation distance by 

W^(i/i,zv2) := / Ty^K,^2)rfPM 
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where 



W^{rii,ri2) = inf I / t^^x - y\) dq{x,y) : g is coupling of r?i , ??2 



denotes the usual L^-Wasserstein distance - w.r.t. the distance '&{\x — y\) - between (not 
necessarily normalized) measures ??i,r/2 £ of equal total mass. 

Lemma 5.5. (i) For any triple of measure-valued random variables v'^v*^^* '■ ^ ~^ M{W^) 
with vfi^'^) = i/^(M'^) = (M'') for a.e. uj ^ Q we have the triangle inequality 

W^(z.i, 1/3) < W^(z.i, 1/2) + W^(z/2, i^s). 

(ii) For each countable families of pairs of measure- valued random variables v'^^v'k • ^ ~^ 
M{W^) with = z^^fe(M'^) for a.e. w G and all k we have 



\ k k J k 



Proof. Gluing lemma (cf. |Dud02| or |Vil09| . chapter 1) plus Minkowski inequality yield (i); (ii) 
is obvious. □ 



For each bounded measurable A (ZW^ let us now define a random measure v\: Q. ^ M.{W^) by 

Note that - by construction - the measures and 1a m'^ have the same total mass. The modified 
transportation cost is defined as 



Ca{!-^) = inf< / c{x,y)dq{x,y) : g is coupling of and 1 a ^'^ ^ = W^{v%Ia^i'^)- 



Put 



with Bn = [0, 2")^^ as usual. 



5.2.2 Semi-Subadditivity of Modified Cost 

The crucial advantage of this modified cost function Q.a is that it is semi-subadditive (i.e. sub- 
additive up to correction terms) on suitable classes of cuboids which we are going to introduce 
now. For n G Nq, /c G {1, . . . , d} and i G {0, 1}'^ put 

:= [0, rf X [0, + 2" • (n, . . . , ifc, 0, . . . , 0) 

These cuboids can be constructed by iterated subdivision of the standard cube B^+i as follows: 
We start with Bn+i = [0, 2""^^)'^ and subdivide it (along the first coordinate) into two disjoint 
congruent pieces B^^l-^ = [0, 2") x [0, 2"+i)'^-i and B^^l^ = B^^l-^ + 2" • (1, 0, . . . , 0). In the A:-th 
step, we subdivide each of the -B^+i = for i G {0, 1}*^^-'^ along the k-th. coordinate 

into two disjoint congruent pieces and ^^*jJ_'_^--'*'=-i'^). After d steps we are done. 

Each of the -B^+i for i G {0, l}*^ is a copy of the standard cube Bn, more precisely, 

B^j^i = Bn -\- 2^ ■ i. 
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Lemma 5.6. Given n G No, A; e audi G {0,1}'= = S^^'f Di = 

^^n...,ifc_i,i) D = Do u L>i = B^^'i'-'^'^-i^ Then 

(z/i)o + i/D, , ud) < 2-("+i)0(2"+i)2'^/2(n+i)-fe/2_ 

Proof. Put Zj(uj) := fi"{Dj) for j G {0,1}. Then Zo,Zi are independent Poisson random 
variables with parameter oq = = ^(Dj) = 2'^("'+^)~'= and Z := /x(Z)) = Zq + Zi is a, Poisson 
random variable with parameter a = 2'^("'+^)~'=+^. 

The measure vd has density ^ on D whereas the measure i>D '■= i^Dq + ^Di has density ^ on 
the part Dq C D and it has density on the remaining part Di G D. If Z = nothing has to 
be transported since P already coincides with ly. Hence, for the sequel we may assume Z > 0. 
Assume that Zq > Z^. Then a total amount of mass ^o-Zi ^ uniformly distributed over Dq, will 
be transported with the map 

T : (xi, . . . , Xk-i,Xk, Xk+i, ...,Xd)<-^{xi,..., Xjk_i, 2**+^ - Xk, Xk+i, ...,Xd) 

from Do to Di. The rest of the mass remains where it is. Hence, the cost of this transport is 

\Zo — Z^ 
2" 

Hence, we get 



^ • 2-" / ??(2'^+^ - 2xfc) dxfc = 2-("+2)g)(2'^+^) • |Zo - Zi|. 



{vb.vd) = 2-("+2)0(2n+i).E[|Zo-Zi|] 

< 2-("+^)6'(2"+i)-E[|Zo-ao|] 

< 2"("+^)0(2"+^) • a^^ = 2-("+i)0(2'*+^)2'^/^("+^)"^/^. 



□ 



Proposition 5.7. For aZ/ n G N and arbitrary dimension d it holds 

Cn+l < c„ + 2'='/2+1.2-(n+l)(d/2+l)0(2"+l). 

Proof. Let us begin with the trivial observations 
and 

wJib„^,/x, ^bA^ E W41b«/x,z^b^)=2'^-W^(1s„/x,i^bJ = 2'^(»+^)-c„. 

\ je{o,i}d y je{o,i}'' 

Hence, by the triangle inequality for an upper estimate for Cn+i — c^i will follow from an 
upper bound for W{f (j2ie{o,i}d ^Bi, , '^B„+i)- 

In order to estimate the cost of transportation from i^^^^ := X^jg^o i}"* ^(o) '■— ^Sn+i 
fixed n G Nq, we introduce {d — 1) further ('intermediate') measures 

ie{o,i}'= 

and estimate the cost of transportation from z/(-fc) to i'(^k-i) for k G {1, . . . , d}. For each k, these 
cost arise from merging 2'="^ ;?airs of cuboids into 2'="^ cuboids of twice the size. More precisely, 
from moving mass within pairs of adjacent cuboids in order to obtain equilibrium in the unified 
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cuboid of twice the size. These cost - for each of the 2^^ ^ pairs involved - have been estimated 
in the previous lemma: 

for A; G {1, ... , d} (and arbitrary i G {0, l}'^"'^). Thus 

d 

k=l 

d 

< ^ 2*^/^ • 2~("+^)0(2"+'^)2'^/^*^"+-'^) 

k=l 

< 4 . 2("+2)('i/2-l) . 0(2'*+^) 

which yields the claim. □ 
Corollary 5.8. // ^ 2-("+i)W2+i)0(2"+i) < oo, we have 

n>l 

?oo := lim Cn 

n— >-oo 

exists and is finite. 

Proof. According to the previous theorem 

lim c„ < CAT + y 2-('"+^)W2+^)0(2"*+^), (5.1) 

m>N 

for each TV G N. As the sum was assumed to converge the claim follows. □ 

5.2.3 Comparison of Costs 

Proposition 5.9. For all d & N and for all n G No 

Cn < c„ + • £(2"). 

Proof Let a box B = = [0, 2")^^ for some fixed n G No be given. We define a measure-valued 
random variable : O ^ A1(M'') by 

with a randomly scaled box B{uj) = [0, Z{cjy/'^)'^ C and Z(a;) = fi'^iB). Recall that Z is a 
Poisson random variable with parameter a = 2"*^. Moreover, note that 

and that \% < & for each a; G Each coupling of A"^ of Is/x"^, therefore, is also a semicoupling 
of 2 and IbI^'^- Hence, 

2"'^-c„<W^(Ab,1b/x). 

On the other hand, obviously, 

2'^'^-c„ = W^(z.b,1bm) 

and thus 

2"''-(c„ - <W^(i/B,AB). 
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If Z > a a transport Tj.i'b = can be constructed as follows: at each point of B the portion ^ 
of vb remains where it is; the rest is transported from B into B\B. The maximal transportation 
distance is \/d ■ Z^^'^. Hence, the cost can be estimated by 

^(Vd-Z^/'^) -{Z-a). 



On the other hand, if Z < a in a similar manner a transport T^Xb = i^b can be constructed 
with cost bounded from above by 

^[^/d-a^''^) ■ (a-Z). 
Therefore, by definition of the function e(.) 



< e \ a 



Z — a\ 
(ZVa)^/^ ■ \Z 



a 



< e(a^/'^yVd-E[Z + a]^/^-E 

1/2 



\Z — a\ 



1/2 



e(2")- Vd' 



2 . . 2""^ 



This finally yields 



c„ < 2 



-nd 



2d. 



□ 



Theorem 5.10. Assume that 



then 



1 r 



l+d/2 



dr < oo 



(5.2) 



Coo — ^oo ^ oo. 



Proof. Since 



v,l+d/2 



dr < OO 



^ 0(2") 

Z.^ 2n(l+fi/2) 
n=l 



< OO, 



Corollary |5.8| applies and yields Coo < oo. Moreover, since '& is increasing, the integrability 



condition (5.2) implies that 



e{r) = sup 



s>r S 



.d/2 







as r — )• oo. Hence, Cqo < Coo by Proposition 5.9 



□ 



The previous Theorem essentially says that Coo < oo if ■!? grows 'slightly' slower than r"'/^. This 
criterion is quite sharp in dimensions 1 and 2. Indeed, according to Theorem 5.2 in these two 
cases we also know that Coo = oo if "i? grows like r'^/^ or faster. 



5.3 Estimates for L^-Cost 

The results of the previous section in particular apply to L^-cost for p < d/2 in d < 2 and to 
-L^-cost for p < 1 in d > 3. A slight modification of these arguments will allow to deduce cost 
estimates for cost for arbitrary p > 1 in the case d >3. 

In this case, the finiteness of Coo will also be covered by the more general results of |HP05j . see 
Theorem 1.3 (i). However, using the idea of modified cost we get reasonably well quantitative 
estimates on Cqo- Throughout this section we assume /3 = 1. 
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5.3.1 Some Moment Estimates for Poisson Random Variables 



For p G R let us denote by \p\ the smallest integer > p. 

Lemma 5.11. For each p G (0,oo) there exist constants Ci (p) , C2 (p) , C3 (p) such that for every 
Poisson random variable Z with parameter a > 1; 

(i) E [ZP] < Clip) ■ «^ e.g. Ci(l) = 1, Ci(2) = 4. 

For general p one may choose Ci{p) = [p]^ or Ci{p) = 2P~^ ■ ([p] — 1)!. 

(ii) E [Z-P ■ l^z>o}] < C2{p) ■ a-P. 

For general p one may choose C2{p) = ([p] + 1)!. 



(Hi) E [{Z - a)P] < CM ■ aPl\ e.g. C^{2) = 1, Ci(4) = 2. 
For general p one may choose C3 = 1p^^ ■ {'2\'^ — 1)!. 

Proof. In all cases, by Holder's inequality it suffices to prove the claim for integer p G N. 



(i) The moment generating function of Z is M{t) := E[e*^] = exp (Q(e* — 1)). For integer p, the 
p-th moment of Z is given by the p-th derivative of M at the point t = 0, i.e. E [ZP] = M^p\o). 
As a function of a, the p-th. derivative of M is a polynomial of order p (with coefficients depending 
on t). As a > 1 we are done. 

To get quantitative estimates for Ci, observe that differentiating M(t) p times yields at most 
2P^^ terms, each of them having a coefficient < {p — 1)1 (if we do not merge terms of the same 
order). Thus, we can take Ci = 2P~^ ■ (p — 1)!. 

Alternatively, we may use the recursive formula 



for the Touchard polynomials T„(a) := E[Z"], see e.g. |Tou56| . Assuming that Tk{a) < [ka)^ 
for all /c = 1 , . . . , n leads to the corresponding estimate for k = n + 1. 

(iii) Put p = 2k with integer k. The moment generating function of {Z — a) is 

Nit) := exp (a(e* -l-t))= exp (f i^M*)) = 1 + f i'M*) + \{'^)W{t) + \{^)h'h'{t) + ... 

with h{t) = ^(e* — 1 — t). Hence, the 2A;-th derivative of N at the point t = is a polynomial 
of order k in a. Since a > 1 by assumption, E[(Z — a)'^^\ = N^'^^\Q) < C^ ■ for some C3. 
To estimate C3, again observe that differentiating N{t) (2k) times yields at most 2^^~^ terms. 
Each of these terms has a coefficient < {2k — 1)! (if we do not merge terms). Hence we can take 
C73(2A:) = 2^^-^ ■ {2k-l)\. 

(ii) The result follows from the inequality 

1 ^ {k + l)\x\ 
~ {k + x)\ 

for positive integers k and x. The inequality is equivalent to 
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For fixed k the latter inequality holds for x = 1. If x increases from x to x + 1 the right hand side 
grows by a factor of {^^)''~^^ and the l.h.s. by a factor of ^'^^'^^ ■ As {x + k + l)x'' < (x + l)^^^, 
the inequality holds. Then, we can estimate 



E 



" 1 






< E 



(fc + 1)! 



(Z + l) 

oo 



E 



a-' 



(Z + k) 
{k + l)l 



(k + l) 



• e 



a" 



E 



a- 



j+k 



\{3 + k)\ 



< 



{k + l)\ 



a" 



If we choose k = \p\ this yields the claim. 



□ 



5.3.2 LP-Cost for p > 1 in d > 3 

Given two measure valued random variables v'^v*:^^ JV{.{W^) with uf{W^) = i/^(M'^) for a.e. 
a; G ri we define their L^-transportation distance by 



Jn 



i/p 



where 



Wp{rji,r]2) = inf 



\x — d6{x, y) 



i/p 



: 9 is coupling of r/i , 172 



denotes the usual L^-Wasserstein distance between (not necessarily normalized) measures r/i, 772 G 
M{W^) of equal total mass. Note that Wp(j^i,z^2) is not the L^-Wasserstein distance between 
the distributions of v' and v'- The latter in general is smaller. Similar to the concave case the 
triangle inequality holds and we define the modified transportation cost as 



inf 



Put 



\x — ?/|^dg(x, y) : q is coupling of z/^ and 1a /^"^ 



2-nd . jg 



-Bn 



with Bn = [0, 2")'^ as usual. 



^(il,...,ifc_i,0) 



Lemma 5.12. Given n G No, A; G {!,..., d} and i G {0,1} put Dq — j^^j^i 

some constant ki depending only on p: 



One may c/ioose = ^2"^ • C3(2p) • C2(2(p - 1)). 



Proof. The proof will be a modification of the proof of Lemma 5.6 An optimal transport map 
T : D ^ D with T^,i'£) = is now given by 



T : {xi, . .., Xk-i,Xk, Xk+i, . . . ,Xd) {xi, . . . , Xk-i, • Xk, Xk+i, ■■■ ,Xd) 



on Do and 



T : (xi, . . . , Xk-i,Xk, Xk+i, ...,Xd)^ (xi, . . . , Xfc-i, 2"+^ - (2"+^ - Xk) ■ — , x^+i, . . . , x^) 

on Di. (If p > 1 this is indeed the only optimal transport map.) The cost of this transport can 
easily be calculated: 

f2" 



'Do 



\T{x) - x\P di?ix) = Zo ■ 2-"" 



2Zo 
Z 



Xk Xk 



dxi 



P+l 



■Zo 



Zo — Z\ 
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and analogously 



\T{x) — x\^ du{x) 



<2np 
P+1 



Zq — Z\ 



z 



Hence, together with the estimates from Lemma 5.11 this yields 

2 rap 



< 



< 



< 



p + l 

P+l 
2(n+l)p 

P + l 

2(™+i)p 



E 



\Zo-Zi\P ^ 

^^Ti 1{Z>0} 



E \Zo-Zi\'i 



■E\\Zo-ao 



2pl 1/2 



E 



^-2(p-l) . ^ 



{Z>0} 



2pl 1/2 _ 



E 



^-2(p-l) . ^ 



1/2 

1/2 



{Z>0} 



P + l 

< Ki ■ 2('^+l)(P+'='-P'^/2) . 2*=(P/2-l)+l 



which is the claim. 



□ 



With the very same proof as before (Proposition 5.7), just insert different results, we get 



Proposition 5.13. For all d G N and all p > 1 there is a constant K2 = K.2ip,d) such that for 
all n G No 



-i/p 



In particular, 



< 71/p 



+ K2 



2-(n+l)(d/2-l) 



1 _ 2-W2-1) ■ 

One may choose K2{p, d) = ki{pY/p ■ Yfk=i '^^''^ < /«i(p)^/^ • 2^^/2+2. 
Corollary 5.14. For all d > 3 and all p > 1 

"Coo := lim < 00. 

n— >oo 

More precisely, 



ci/p 



1/p AKi{py/p 

- '0 +2-i-2-<i/2- 



Comparison of costs "c^ and Cn now yields 

Proposition 5.15. For all d> 3 and all p > I there is a constant K3 such that for all n G No 

Cl/P < ?l/p + ^3.2"(l-'^/2). 



Proof. It is a modification of the proof of Proposition 5.9 This time, the map T : B ^ B 



T : 



a 



defines an optimal transport T*z^b = ^B- Put t' = T'{d,p) = ^-^^ \x\'p dx. (This can easily be 
estimated, e.g. by t' < j^c^^^^ if P > 2.) The cost of the transport T is 



|T(x) — x\P du^ix 



t' • 2"p • Z 



< t'-2"P-Z- 



Z 



a 



i/d 



Z 



a 
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The inequality in the above estimation follows from the fact that \t—l\ < \t—l\-{t +. . .+t+l\ = 

jt*^ — 1| for each real t > 0. The previous cost estimates holds true for each fixed u (which for 
simplicity we had suppressed in the notation). Integrating w.r.t. dF{uj) yields 





z 


p- 


z- 


— 1 






a 





< t' ■ 2''P ■ a-P ■ a ■ Cs ■ aP/^ = 

A/p. _ ^l/p* < ^3.2n(l-<i/2). 



and thus 



a 



2pl 1/2 



^p . 2"(a!+P-'ip/2) 



□ 



Corollary 5.16. For all d>3 and all p > 1 



Coo — 



< OO. 



5.3.3 Quantitative Estimates 

Throughout this section, we assume that ??(r) = with p < p{d) where 



p < p{d) :-- 



GO, for d > 3 
1, for d = 2 
\, for d = l. 



Proposition 5.17. Put T{p,d) = • (r(| + l)V<i . tt'^/'^Y . Then 

Coo > Co > T(p,d). 

Proof. The number r as defined above is the minimal cost of a semicoupling between £ and a 
single Dirac mass, say Sq. Indeed, this Dirac mass will be transported onto the d-dimensional 
ball Kr = {x e '^'^ : \x\ < r} of unit volume, i.e. with radius r chosen s.t. S.{Kr) = 1. The cost 
of this transport is \x\p dx = j^r^ = r. 

For each integer Z > 2, the minimal cost of a semicoupling between £ and a sum of Z Dirac 
masses will he > Z ■ r. Hence, if Z is Poisson distributed with parameter 1 



Co > E[Z] ■T = T. 

Remark 5.18. Explicit calculations yield 



□ 



r{p, 1) 



1 



l+P 



r{p,2) 



2 + p 



-p/2 



r{p,S) 



3 + p 



-) 



p/3 



whereas Stirling 's formula yields a uniform lower bound, valid for all d & N ( which indeed is a 
quite good approximation for large d ) 



r{p,d) > 



d ( ^y'^ 



d + p \27re / 

Proposition 5.19. Put 9 = T{d,p) = /jq /jq \x — y\P dydx. Then 



< Co < f. 



Moreover, t < (x_^p)(\_)_p/2) ' ^^^^ f^'^ allp>2 and f < (g) for all < p < 2 
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Proof. If there is exactly one Poisson particle in Bq = [0, l)'^ - which then is uniformly distributed- 
then the transportation cost are exactly T{d,p). In general, r still is an upper bound for the 
cost per particle. The number of particles will be Poisson distributed with parameter 1. The 
lower estimate for the cost follows from the fact that with probability there is exactly one 
Poisson particle in Bq = [0, l)*^. 

Using the inequality (xf + . . . + x^)^/^ < d^^"^'^ ■ (x^ + . . . + x^) ~ valid for all p > 2 - the upper 
estimate for r can be derived as follows 



/ / \x-y\Pdydx < dP/^-iy^ f f 



Jo Jo 

1 



[0,1]'* 

s - t\Pdsdt 
■d^l\ 



\xi — Uil^ dydx 



(l+p)(l+p/2) 

Applying Holder's inequality to the inequality for p = 2 yields the claim for all p < 2. 
Theorem 5.20. For all p < 1 and d > 2p 



□ 



d 



p/2 



d + p \27re J 
whereas for all p > 1 and d >3 

d V/^ 

' \2TTe J 



< 



1 





ey + (p + i)(2"!/2-P_i) 



d + p 



< clip < 



dV2 



61/2 A[(l+p)(l+p/2)]i/p 



i7^ + 28--i 



i/p 



Proof. Proposition 5.17 and the subsequent remark yield the lower bound 

d ( ^Y'^ 



d + p \ 27re J 



< T < Cn 



valid for all d and p. In the case p > 1 the upper bound follows from Proposition 5.19 and 

I/P ^1/2 



Corollary 5.14 by 



+ 



4k 



< 



2-d/2 - 61/2 ^ [(1 + + p/2)]i/p 



i7^ + 28-«i 



i/p 



In the case p <1, estimate (5.1) with &{r) 



p+1 



Coo < Co + ^ 2 



{m+l)(d/2+l) 



1 



,(m+l)(p+l) 



Co + 



m=0 



p+1 ' (p+l)(2rf/2-P-l) 

provided p < d/2. Together with Proposition 5.9 this yields the claim. 
Corollary 5.21. (i) For all p € (0, oo) 

1 cU." cj/^ 1 

< lim inf < lim sup —r^ < 



□ 



/27re 



dV2 



V6A[(1+p)(1+p/2)]Vp' 



Note that the ratio of right and left hand sides is less than 5, - and for p < 2 even less than 2. 
(a) For all p G (0,oo) there exist constants k,k' such that for all d > 2(p A 1) 

k-dP/^ < c^ < k'-dP/^. 
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6 Appendix. Optimal Semicouplings with Bounded Second Marginals 



The goal of this chapter is to prove Theorem 2.1 (= Theorem 6.6), the crucial existence and 
uniqueness result for optimal semicouplings between the Lebesgue measure and the point process 
restricted to a bounded set. The theory of optimal semicouplings is a concept of independent 
interest. Optimal semicouplings are solutions of a twofold optimization problem: the optimal 
choice of a density p < 1 oi the first marginal fii and subsequently the optimal choice of a 
coupling between pfii and fi2- This twofold optimization problem can also be interpreted as a 
transport problem with free boundary values. 

Throughout this chapter, we fix the cost function c(x, y) = i^dx — y\) with ~ as before - being 
a strictly increasing, continuous function from to M+ with ??(0) = and lim i?(r) = oo. 

r— >oo 

Lemma 6.1. Given a finite set H = {^i, . . . , ^k} C M'* and a probability density p G L^(K'^, £,). 
(i) There exists a unique coupling q of p£ and ^ = ^ X^ge= '^C which minimizes the cost function 
Cost(-). 

(a) There exists a (Sj-a.e. unique) map T : {p > 0} — J- H with T^,{p2,) = a which minimizes 
f c{x,T{x))p{x)d£{x). 

(Hi) There exists a (2,-a.e. unique) map T : {p > 0} — )• H with T^{p£) = a which is c-monotone 
(in the sense that the closure of {{x,T{x)) : p{x) > 0} is a c-cyclically monotone set), 
(iv) The minimizers in (i), (ii) and (Hi) are related by q = {Id,T)^,{pS.) or, in other words, 

dq{x,y) = ddT(x)iy) P{x)d£{x). 

Proof. We prove the lemma in three steps. 

(a) By compactness of Il{pil,cr) w.r.t. weak convergence and continuity of c(-, •) there is a 
coupling q minimizing the cost function Cost(-) (see also |Vil09j . Theorem 4.1). 

(b) Write p£ =: A = Yli=i ^« where Aj(.) := q{. x {Ci}) for each i = 1, . . . , k. We claim that the 
measures (Aj)j are mutually singular. Assuming that there is a Borel set N such that for some 
i ^ j we have Xi{N) = a > and Xj{N) = /3 > we will redistribute the mass on N being 
transported to and in a cheaper way. This will show that the measures (Aj)j are mutually 
singular. In particular, the proof implies the existence of a measurable c-monotone map T such 
that q = {Id,T)^{p£). 

W.l.o.g. we may assume that {p£,){N) = a + (3. Otherwise write p = pi + P2 such that on N 
dqi{x) + dqj{x) = d(/>i£)(x) and just work with the density pi. 

Put f{x) := c{x,^i) — c{x,S^j). As c(-,-) is continuous, f is continuous. The function c{x,y) 
is a strictly increasing function of the distance |j; — y|. Thus, the level sets {/ = 6} define 
(locally) (d — 1) dimensional submanifolds (e.g. use implicit function theorem for non smooth 
functions, see Corollary 10.52 in [Vil09j ) changing continuously with b. Choose 6o such that 
p£({/ <bo}nN) = a (which implies p£.{{f > bo} n N) = (3) and set Ni := {/ < 60} n iV and 
Nj := {/ > 60} n N. 
For I = i,j 

d\i{x) := d\i{x) - lN{x)dXi{x) + lN^{x)d{p£,){x). 

For I i,j set A; = A;. By construction, q = X]f=i ^ ^£,1 is a coupling of p2 and a. Moreover, 
q is c-cyclically monotone on N, that is Vxi G Ni,Xj G Nj we have 

c{Xi,^i) +c{Xj,^j) < c{Xj,^i) +c{Xi,^j). 

Furthermore, the set where equality holds is a null set because c{x, y) is a strictly increasing 
function of the distance. Then, we have 

Cost(g) — Cost(g) 

c{x,(,i)dqi{x) + c{x,^j)dqj{x) - / c{x,S,i)dqi{x) - c{x,^j)dqj{x) > 0, 

N JNi JNj 
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by cyclical monotonicity. This proves that Aj and Xj are singular to each other. 
Hence, the family {\i)i-i„,,^k is mutually singular which in turn implies that there exist Borel 
sets Si C M'^ with (j-Si = R'^ and Xi{Sj) = for ah i ^ j. Define the map T : M'^ ^ H by 
T{x) := for all x E Si. Then q = {Id,T)4p£). 

(c) Assume there are two minimizers of the cost function Cost, say qi and q2- Then ^3 := 
2(91 + 92) is a minimizer as well. By step (b) we have qi = {Id,Ti)^p2, for i = 1,2,3. This 
implies 

d5T3{x){y) dp2,{x) = dq3{x,y) = d (^^qi{x , y) + ^q2ix , y) 

= d (^^5T^^^){y) + ^5T^(^^'){y)^ dp£{x) 

This, however, implies Ti(x) = ?2(x) for p£ a.e. x G M"* and thus qi = q2- D 
Remark 6.2. In the case ^{r) = r^, there exists a convex function : {p > 0} — t- M such that 

T{x) = Vif{x) for 2,-a.e.x. 

2-p 

More generally, if'd{r) = r^ withp > 1 then the map T is given as T{x) = x+\\/ip{x)\p~''^ ■Vip{x) 
for some -convex function : {p > 0} — ?■ M. 

Proposition 6.3. For each finite set H C M'^ there exists a unique semicoupling q of £ and 
a = X^ggs ^5 which minimizes the cost functional Cost(-). 

Proof, (i) The functional Cost(.) on A4(M.'^ x M'^) is lower semicontinuous w.r.t. weak topology. 
Indeed, if — )• weakly then with Cfc(x, y) := min{i?(|x — y|), k} 



liminf Cost(ry„) > sup 

" k 



n 



lim / Ckdr]r, 



sup / Ck dr] = Cost{r]) . 

k J 



(ii) Let Q. denote the set of all semicouplings of £ and a and £2i the subset of those g G £} 
which satisfy ^Cost{q) < inf^/gQ Cost(g') =: c. Then Qi is relatively compact w.r.t. the weak 
topology. Indeed, q{]R.'^ x CH) = for all (7 G 0i and 

q{CKr{E) X H) < . Cost(g) < 

for each r > where K^C^) denotes the closed r-neighborhood of H in M'^. Thus for any e > 
there exists a compact set K = Kj-i^) X H in M'' X M'^ such that g(C-ftr) < e uniformly in g G £2i. 

(iii) The set £} is closed w.r.t. weak convergence. Indeed, ii qn ^ q then (7ri)^,(7„ — )• (vri)^,g and 
{'^2)*qn {T^2)*q- 

Thus, 111 is compact and Cost(.) attains its minimum on (or equivalently on lli). 

(iv) Now let a minimizer q of Cost(.) on £l be given and let A = (7ri)*(/ denote its first marginal. 
Then A = p • £ for some density < p < 1 on M'^. Our first claim will be that p only attains 
values and 1. 

Indeed, put U = {p > 0}. According to the previous lemma 6.1, there exists an a.e. unique 
'transport map' T : U ^ 3 s.t. 

q={Id,T%X. 

For a given 'target point' ^ G H, C/^ := [/ Pi T^^{^) is the set of points which under the map 
T will be transported to the point ^. Within this set, the density p has values between and 
1 and its integral is 1. If the density is not already equal to 1 we can replace it by another 
one which gives maximal mass to the points which are closest to the target ^. Indeed, put 
r(0 := inf{r > : £(Kr-(0 n U^) > 1} and A := p • £ with 
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Then 

q:= {Id,T)^\ 

defines a semicoupling of 2, and a with Cost(g) < Cost(g). Moreover, Cost(q) = Cost(g) if and 
only if p = p a.e. on M'^. The latter is equivalent to p £ {0, 1} a.e. 

(v) Assume there are two optimal semicouplings qi and q2 whose first marginals have density 
Ijji and l[/2, resp. Then q := ^{qi + q2) is optimal as well and its first marginal has density 
+ lu2)- By the previous part (iv) of this proof the density can attain only values or 1. 
Therefore, we have Ui = U2 (up to measure zero sets) and qi = q2- □ 

Lemma 6.4. Given a bounded Borel set A C M**, let Mcount{-A) = {a G Mcount(W^) '■ <y(W^ \ 
A) = 0} denote the set of finite counting measures which are concentrated on A. Define T : 
■M.count{A) — )• A^(M"'xM'^) the map which assigns to each a G M-counti-^) the unique q G ns(£, a) 
which minimizes the cost functional Cost(.). Then T is continuous (w.r.t. weak convergence on 
the respective spaces). 

Proof, (i) Take a sequence ((T„)„ C M.count{A) converging weakly to some a G Mcounti^)- Put 
qn ■= T((Jn) for n G N and q = T{a). Wc have to prove that g„ — >■ q. 

(ii) The weak convergence o"„ — t- cr implies that finally all the measures o"„ have the same total 
mass as a, say k. Hence, for each sufficiently large n G N there exist points x", ■ ■ ■ and Borel 
sets Si,. . . such that 

k k 

an = ''^Sxf, qn = Igf ^ ^xf- 
i=l i=l 

Similarly a = Yli=i and q = Yli=i ^Si^ ® ^xi with suitable points xi, . . . , and Borel sets 
Si,.. . ,Sk. Weak convergence moreover implies that for each i = 1, . . . ,k 

— )■ as n — >■ 00. 

(iii) Based on the representations of q and (T„, we can construct a semicoupling g„ of Q and (7„ 
as follows 

k 
i=l 

Then by continuity of •& and dominated convergence theorem 

lim sup Cost(g„) = lim sup V] / '^ilv — x2\)dy = / i^Hy — Xi\)dy = Cost{q). 

And of course Cost(g„) < Cost(g„). Thus 

limsupCost((;„) < Cost(g). 

n 

(iv) The sequence {qn)n is relatively compact in the weak topology of M.{W^ x W^). Therefore, 
there is a subsequence, denoted again by {qn)n, converging weakly to some measure q G A4{W^ x 
M*^). It follows that {'K2)*(ln {'^2)*Q and thus (712)*^ = a. Similarly, (7ri)*g < £. Thus 
q G Ils{£,a). Lower semicontinuity of the cost functional implies 

Cost(g) < liminf Cost(q„). 

n->oo 

(v) Summarizing, we have proven that g is a semicoupling of ii and a with 

Cost(g) < Cost(g). 
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Since q is the unique minimizer of the cost functional among all these semicouplings, it follows 
that q = q. In other words, 

lim T((T„) = T( lim c7„). 

n— >-oo n— >-oo 

This proves the continuity of T. □ 



For a given u let us apply the previous results to the measure 

cr = Iah'^ = ^ 

for a realization ji'^ of the point process. Then, there is a unique minimizer - in the sequel 
denoted by q^ - of the cost functional Cost among all semicouplings of 2, and lAfJ-^- 

Lemma 6.5. For each bounded Borel set ^ C M'^ the map u q"^ is measurable. 



Proof. We saw that the map T : Mcount{A) — )• M{R x M"), a i-> T(a) assigning to each 
counting measure a its unique minimizer of Cost(.) is continuous. By definition of the point 
process, co /j!^ is measurable. Hence, the map 

is measurable. □ 

Theorem 6.6. (i) For each bounded Borel set A cW^ there exists a unique semicoupling Qa 
of 2, and (1a/x*)P which minimizes the mean cost functional £ost(.). 

(ii) The measure Qa can be disintegrated as dQA{x,y,uj) := dq^{x,y) dF{uj) where for ¥ -a. e. oj 
the measure q^ is the unique minimizer of the cost functional Cost(.) among the semicouplings 
of Q and IaA*'^- 

(Hi) (t05t{QA) = J^Cost{q'X)dF{oj). 

Proof. The existence of a minimizer is proven along the same lines as in the previous proposition: 
We choose an approximating sequences Qn in M{^^ xW^ x Q,) - instead of a sequence g„ in 
Ai{W^ X W^) ~ minimizing the lower semicontinuous functional €o5t{.). Existence of a limit 
follows as before from tightness of the set of all semicouplings Q with (tost{Q) < 2infg (Lo5i{Q). 
For each semicoupling Q of ii and /x*P with disintegration as g*P we obviously have 



(tost{Q)= / Cost{q'^)dF{oj). 
Jn 



Hence, Q is a minimizer of the functional £ost(.) (among all semicouplings of £ and /x*P) if 
and only if for P-a.e. a; G the measure q'^ is a minimizer of the functional Cost(.) (among all 
semicouplings of ii and /x'^). 

Uniqueness of the minimizer of Cost(.) therefore implies uniqueness of the minimizer of £ost(.). 

□ 

Corollary 6.7. For each z € M*^ and each bounded Borel set A C the push forward of the 
measure Qa under the translation Tz ■ {x,y,u) ^ {x + z,y + z,u) ofW^xW^xU, coincides 
with the measure Qz+A- 

Proof. Since ii is invariant under the translation x ^ x + z and /x*P is invariant under the 
translation (y, uj) ^ [y + z, oj) the claim follows from the uniqueness of the minimizer of the 
cost functional £ost(.). □ 
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Remark 6.8. As before for a finite set H C M*^ put a = X^^gst^^- Let q be a semicoupling 
of £ and a. Then, q minimizes Cost(.) iff the support of q is c-cyclically monotone and q is 
c- sequentially monotone in the following sense: 

n n 

i=l i=l 
for all n G N, {{xi,^i)}f^^ £ supp(g), Vx„+i ^ supp((7ri)*g). 

Proof. Let q be the Tinique minimizing semicoupling. The cyclical monotonicity follows from 
the general theory of optimal transportation. Put U := supp((7ri)*g). Assume that q is not 
sequentially monotone. Then, there are n G N,a; = Xn+i € Ct/, {{xi,^i)}f^i G supp(g) such that 

n n 
^c(Xi,^i) > ^c(Xj+l,^i). 

i=l 1=1 

By continuity of the cost function, there are (compact) neighborhoods Ui of Xi and Vi of such 
that Un+i n [/ = and 

n n 
i=l i=l 

whenever Ui G Ui and vj G Vj. Moreover, as supp((T) is discrete we can assume (by shrinking 
Vj slightly if necessary) that Vj fl supp(cr) = {^j}. As {xi,^i) G supp(q') for 1 < z < n we 
have inUqiUi x fe}) > 0. Set A := M{q{Ui x {^i}) , . . . , q{Un x {^n}) , &{Un+i)} ■ Then, we 
can reallocate mass to define a new measure with less cost. Indeed, we can choose subsets 
Ui C Ui, Ui X C supp(g) with £y{Ui) = A and define a new measure q by 

^ n \ ^ 

dq{x, y) = dq{x, y) " " ^ ^C/ixUjl^' y)d^{x) + " 1^ lc/,+ixte}(^' y)d£(x). 

1=1 i=l 

By assumption, we have Cost(g') < Cost(g). Hence, q is not minimizing Cost. 

For the other direction let us assume that q is cyclically monotone and sequentially monotone 
but not minimizing Cost(.). Then, there is a Borel set U ^ U{= supp((7ri )*</)) (by uniqueness 
of optimal transportation of fixed measures) and a unique Cost minimizing coupling q of l^Q 
and a such that Cost(g) < Cost(g) and the support of q is cyclically monotone. As U ^ U there 
is some z G U\U which is transported by q to say. For ^ £ S set 5^ := {x G M'^ : {x,^^) G 
supp(g)} and similarly for q. By sequential monotonicity of q for all xq G S^^ we must have 
c(xo,Co) < c(z,^o)- Moreover, the set {x G S^^ : c(a;,^o) = c{z,S,o)} is a £ null set. Thus, there 
is a set S^g C S^^ of Lebesgue measure one such that for all x G S^^ we have c(x,^o) < c{z,Co)- 
By the first part, we know that a minimizing semicoupling is sequentially monotone. Thus, 
S^Q C U and also S^^ C U (in particular ii S = {^o} we are done). 

Moreover, by assumption there is some xi G S^q\S^q which is transported by q to some G H. 
Then, S^^\S^^ is not empty. If S^^ fl CC/ 7^ we choose X2 G S^-^ n CC/ and stop. If S^^ C U there 
is X2 G S'^^\S'^j which is transported by q to some ^2- If ^2 G {Coi Ci} (that is ^2 = ^0) wc choose 
X2 G 5*^2 ^ ^(1 stop. Otherwise we proceed in the same manner until cither n CU 7^ 
or ^fc G {^0) • • • )^fc-2}- By this procedure, we construct a sequence xq, . . . ,Xk such that Xj G 
S^.nS^._^ ioT 1 < j < k-l, xo e S^q\U and either Xk G %\C/ or Xk G % n52/^_^ = S^.nSyi^_^ 
for some < j < A; — 2. In the latter case, we have by cyclical monotonicity for q and q 

k k k 

^c{xi,^i) < Y^c{Xi+l,^i) < ^c{Xi,Ci), 
i=j i=j i=j 
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where = and Xk+i = Xj. Hence, we have equality everywhere. However, we can move the 
Xi shghtly to get a contradiction. Thus, we need to have x^ G S^^\U. Then we have by the 
sequential monotonicity of q and q 

k-l k-l k-l 

^c{Xi,^i) < ^c(Xi+l,^j) < '^c{xi,Ci). 
1=0 2=0 i=0 

Hence, we need to have equality and therefore a contradiction as before. Hence, q = q. □ 
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